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Abstract: The Luttinger model was introduced to illustrate the theory of Tomon- 
aga via an exactly soluble model. It became soon the subject of great interest also 
on the part of Mathematical Physics and a key to the investigations of the mathe- 
matical properties of Condensed Matter Physics. This paper reviews aspects of the 
above developments relevant for renormalization group methods. 

§1. The Luttinger model. 

The model describes many body systems ol two kinds of fermions on a line whose 
fields are V'S x j with a = ± specifying the creation and annihilation operators 
for fermions located at a point x S [—L/2, L/2] C and distinguished by the 
label a; = + or a; = — . The Hamiltonian, "(kinetic energy) + (chemical poten- 
tial) + (potential energy)", is written as 



(1.1) 

+ \ I dxdyTjjl .^'ip^ .^-tp+ y'ip_^yV{x - y) 

J -L/2 

where A is a coupling constant and v{x — y) is a smooth short range pair potential 
{e.g. v{x — y) = Q \x — y\ > Pq^ for some "range" Pq^) . The parameter vp is 
the "velocity at the Fermi surface", vp = pp/m if m is the mass of the particles: 
we have set h = m = 1. 

This model was introduced in [Lu63] and it illustrated Tomonaga's theory of spin 
fermions in one dimension, [To50], which shows the remarkable phenomenon 
of the anomaly of the ground state: i.e. a ground state with a density of states 
which does not have a discontinuity at the Fermi momentum k = pp but its 
graph becomes vertical with a vertical tangency exponent a(A) — O(A^) called the 
anomaly" of the Fermi surface. 

The model supposes a priori that the "physically significant" part of the Hamil- 
tonian is described by quasi particles. This means realizing that the Schwinger 
functions of the ground state of a spin Fermi gas can be identified, to leading 
non trivial order in the coupling constant A, by thinking that the system in fact 
consists of two particles with energy close to the Fermi energy Pp/2, one with mo- 
mentum close to -\-pF and the other with momentum close to —pp whose fields, 
in a suitable superposition, yield the field tp^ of the observable particles. 

Writing a momentum close to ipp as ±pp + k the free field is expressed as 

^ E<-^^'"""' (1-2) 

LO = ± 

where calling af^ (^t+ppuj uji if Qfp the creation and annihilation operators 
for the Fermi particles, the quasi particles fields are 

<. = ^Ee^''^<fc (1-3) 
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Hence if the operator Tq 



OJ k>0 



approximates the free fcrmionic kinetic energy plus the chemical potential energy 
for fc ~ 0, then the quasi particle fields at time t will be ipux u ~ e^*-^°*V'^^e'F*^o* 
= it with 



- e*^°^^.e-*^° = _L^e±(^'=-+-'=)at (1-5) 



I do not repeat here the heuristic analysis showing that taking the two kinds of 
particles Hamiltonian (1.1) and approximating it with the free kinetic energy (1.4) 
and expressing the fcrmionic field via (1.2) means 

(1) thinking the particles with positive momentum as distinct from the ones with 
negative momentum and, furthermore, 

(2) replacing the dispersion relation e{p) = with p = ±pF + k with e{±, k) = 
p^/2± kpF and, also, 

(3) allowing k to take all values rather than k > —pp or k < pp- 

Luttinger realized that Tomonaga's theory amounted, to leading order in A, to 

(i) replacing the real fermionic particle fields with the ip'^ fields defined by th r.h.s. 
of (1.2) and (1.3), 

(ii) replacing the "real" fermionic particle kinetic energy defined as 



/^/"^ - 1 ~ ■t-^ 1 

^^^dxi,t ^{dl-pDi^- s Y.2^P'-PF)'^t'^P (1-6) 



where p = 2im/L, n = 0,±1,±2, by the operator To above, which can be 
written 



fL/2 

To= dxY,^+J-ivFUjd,)i;-, (1.7) 

J-L/2 -Tr 



and 

(iii) replacing the potential energy by the expression in (1.1), rather than consid- 
ering the usual pair potential which would be written as 

X dxdyi^^i^~ip:^il)~v{x-y) (1.8) 

J-L/2 

Note that comparing (1.8) with the potential energy in (1.1) and with (1.2) several 
"cross terms" involving fields with different quasi particles labels are absent. 

The Luttinger approximation was stated to be exactly soluble in the remarkable 
paper [Lu63] and the solution did yield the expected anomaly of the Fermi surface 
mentioned above, thereby providing a simple explanation of the phenomenon. 

The exact solution was, however, not really correct because of an error on the 
Fock space; canonical commutation relations which, in infinitely many degrees of 
freedom systems, do not have a unique representation: a matter that is now well 
understood but that was not so clear at the time. Nevertheless the attempt in 
[Lu63] contained the important ideas indicating strongly that the model could 
probably be really solved exactly. The exact solution was discovered a little later 
by Mattis and Lieb, [ML65]. The real value of the anomaly and Luttingcr's agree 
to leading order in A {i.e. to second order) but they differ in the higher orders 
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(which vanish identically in Luttingcr's expression while in [ML65] the anomaly 
is analytic in A near A = but it is not a polynomial of second order); see also 
[BGM92]. 

The work [ML65] should be seen as a part of the series of exact solutions of simple 
models which, starting with the 2-dimensional Ising model led to the Bethe ansatz 
(which is not necessary in the case of the Luttinger model) and to exactly solvable 

6 and 8 vertex models, [ML66], [Ba82]. 

The Luttinger model is closely related to the Thirring model of quantmri field 
theory, [Lu63] : for an analysis of the relation between the two models sec [Ma94] . 

§2. Schwinger functions and ground state. 

More recently there has been renewed interest, particularly after the discovery 
of high temperature superconductivity, in understanding the properties of con- 
densed matter from a fundamental point of view, i.e. without relying on heuristic 
arguments, whenever possible. One of the directions in which research has been 
stimulated is the application of multiscale analysis, i.e. of the renormalization 
group methods, to investigate the large distance properties of the correlations in 
the ground states of the simplest systems, i.e. Bose or Fermi gases with weak 
short range interactions: [BG90], [FT90], [BGPS94], [Sh92], [BG95]. Other meth- 
ods have also been applied leading to complementary results valid also at strong 
coupling, for instance [KL73] and [KLY88], [LYOO], or for exact solutions [LL63], 
[LW68]. 

A mathematically rigorous theory of the ground state of the 1-dimensional Fermi 
gas turned out surprisingly hard even at small A, probably because the basic for- 
malism, even today, is still waiting to acquire an established shape. Furthermore 
the case of spinning electrons with attractive rotationally symmetric pair interac- 
tion still defies research. The situation is worse in higher dimension where only 
formal developments seem to be available (see [BG90], [FT90], [Sh92]): there- 
fore higher dimension will be left out of the present discussion, but for sporadic 
comments. 

In dimension 1 the Hamiltonian of a system of spin 0, mass m = 1, fermionic 
particles in a periodic box [— L/2, L/2] (and in the grand canonical ensemble) is 

^ / 1 \ 
H = Y,(-^^^.~f^]+'2>^J2vixi-Xj) (2.1) 

s=l ^ ^ i<j 

with A'^=number of particles or equivalently, in second quantized form, 

f.L/2 
-L/2 



f ' ~ 1 

i^-^) 

+ X dxdy v{x - y) 'ijj+'ijj^ '4>+i>, 

J-L/2 



-L/2 

where ip is the usual Fermi field (rather than its approximation in (1.2), (1.3)). 
If p is the density of the gas the simplest question that one can ask is about 

the behavior near p = ±7rp = ± pp = ±\/2/I of the Fourier transform of the 

one particle reduced density matrix. Note that if A = it is known that the one 
particle reduced density matrix at inverse temperature P = +oo is 

with Fourier transform Xpf{p) = 1 if \p\ < PF and zero otherwise. 
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In the following it will be more convenient to study the more general Schwinger 
functions 



^(71 ■■■(7rL (j^l^l: • ■ • : ^n'^n) 



(-1)'^ lim lim • 

/3— J-OO A^OO 



(2.4) 



Tr e-0H 



where ai = ±1, tt is the permutation of (1, ... , n), such that ^^(i) > ^7r(2) > • • • > 
and (—1)'^ is the permutation parity; n is even. In this way Xpf{p) is the 
Fourier transform of S |-((0~,a;), (0,0)). 

§3. Functional integration for 1— dim spinless fermions. Luttinger— 
Kohn— Ward determination of the chemical potential. 

Denoting ^ — {xo,x),r] — {yo,y) & B?,k= {ko,k) € the Schwinger functions 
introduced in §2 can be usefully expressed as functional integrals 

Sai-an (6, • • • , ^n) = jim Um • 

/3— >oo A— >cx) 

/P(#) I-L%Io 'l't'l'i'l'>v^(^o-yo>(^-y)'iidr, ^ai ...^|„ (3.1) 

where the arc Grassmanian variables", see for instance sec. 3.1 in [BG95]. 
The "integration" with respect to P{dip) is defined on monomials . . . i/i^ by 

assigning the value of (V'J . . . V'^^^) J P{dip) • • • V"^^ via Wick's rule with 
propagators {tp7'ip^) — {'fpt'fp^) — and 



(V;^-^+) = 3(e-o 

9i0= E i-ir9{xo+no,x + nL,p) , 



(3.2) 



5(0 



(27r)<i+i / -i 



^-ipoXQ—ipx 



dpodp 



so that c/(^) is periodic in x with period L and antiperiodic in xq with period /3. 
The integration is extended linearly to even polynomials and to functions of the 
Grassmanian fields that admit an entire even power series expansion, like the 
exponential in (3.1). 

In this way the integrals in the numerator and in the denominator of (3.1) are 
defined as formal power series in A. The series can be quite easily shown to be 
convergent for |A| < 0{L~^) and one of the goals is to find conditions under which 
they can be analytically continued to values of A which are, possibly, small but 
independent of the size L of the system. 

If we insist in fixing a priori the chemical potential the first effect of the interaction 
will be that the singularity of the Fourier transform of {ip^ipQ), see (1.3), will no 
longer be where it is located when A = 0, namely at fc = ±pf = but 
(reasonably) it will be shifted by 0(A). 

It is, therefore, more convenient to keep the location of the singularity fixed at 
a prefixed value pp- This can be achieved by taking a A-dependent chemical 
potential fi = \p^p + v, see (2.1), where v has to be conveniently chosen as a 
function of A. 
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This moans that instead of (3.1) one studies the same expression with the argu- 
ment of the exponential modified into 



-V j I ip+ip,d^-X / ip+ip,ip+ip:^6{xo-yo)v{x-y)d^dri 

J-L/2J0 J-L/2J0 

(3.3) 

In a series of basic papers, [Lu60], [KL60], [LW60], Luttinger, Kohn and Ward 
point out that this determines v as a power series in A and it has the important 
effect of generating a power series for the Schwinger functions which is finite to 
all orders, uniformly in the size of the system: had we fixed the chemical potential 
fi rather than the Fermi momentum pp we would have obtained a power series in 
A with coefficients diverging to all non trivial orders as L ^ 00. 

Of course the latter divergence does not mean that the theory with a fixed chem- 
ical potential cannot be defined: it simply means that such a theory will have 
Schwinger functions with a singularity at a Fermi momentum which is different 
from ^/2i2 so that the expansion with reference to a free field with singularity at 
^/2fl contains diverging expressions. It makes also clear that it is likely to be more 
convenient to develop a perturbation theory of the ground state at fixed Fermi 
momentum rather than at fixed chemical potential. Should one wish to study the 
problem at fixed chemical potential fi, after developing the theory at fixed pp and 
obtaining the chemical potential correction u in terms of X,pf, one can imagine, 
that ijl = pf + fiPF, A) and solve this relation for pp, as a function of /z, A 

A related important result due to Luttinger ( "Luttinger' s theorem") states that 
fixing the Fermi momentum is equivalent to fixing the density which is p = irpp 
whether A = or not: this was shown by Luttinger, [Lu60], to hold to all orders 
of perturbation theory. A mathematically rigorous proof would be desirable. 

§4. The ultraviolet problem. 

If is the range of the interaction potential and e{p) = (p^ —p\)/2 the prop- 
agator (3.2) can be written as sum of two terms 



9(f) = !/'*°'«)+S<^°'«) 

{2'kY J -tqo + e{q) 



(4.1) 



where ^ = {xo,x), n — {qa,q). The term 3^^°^ is the "ultraviolet component of 

the propagator" and g*--*^' is the "infrared component" . 

The decomposition can be used to introduce two auxiliary "independent" Grass- 

manian fields Vj'*'^' and Vg"*^' so that V'j = + This means that in 

evaluating the integrals in (3.1) we can replace V'x with V'^'''^^ -|- and perform 
the Grassmanian integration following the Wick's rule with propagators 

(Vf = - ^) a = (>0),(<0) (4.2) 

while all the other itors vanish. Calling ^(V') the expression (3.3) this is 

also written as the identity ("Fubini's theorem" for Grassmanian integrals) 
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Remark: attention should be paid to the fact that (4.1) deals with the "infinite 

volmnc'' (/? 00, L 00) limit g(^) rather than with the g{£^) of (3.2). Hence one 
should really deal with g and decompose the latter into the sum g^^°''(0+5^~'^H0- 
This generates a great variety of "small" problems both in the ultraviolet and, 
later, in the infrared analysis. There is no few words way out of this (well known) 
difficulty. Here we choose to ignore it except for a few necessary comments when 
needed, because it is discussed widely in the literature, see [BG95], [BMOOa]. 

The idea is to perform first the integral over the "high frequency part" ^p^^'^^ of 
the field both in the mmicrator and in the denominator of (4.3). Fixing attention 
on the denominator (simpler than the numerator which, however, can be treated 
in the same way) the result will be written 

J P(#) e^W = J P(#(^o)) e^'"'^'^*-"') (4.4) 

which in fact is a definition of V^^\ the "effective potential on scale po" ■ 
Technically there is a lot of work behind the latter relation: one has to show that 
the result of the integration of the field tp^''^'> can be written in the form of an 
exponential of an effective potential F'^"^: this means showing that the result can 
be expressed as an exponential of 

m>0,p<m 

and the kernels Vm°p(^i • • • ^m) are 

(i) analytic in A, v. 

(ii) decay exponentially on scale po: which means that they are bounded by 
C„ exp — cPq ^d(^i, . . . ^2n) where d{^i, . . .^2n) is the length pf the shortest path 
connecting the points Ci, • ■ • ^2n in P^" and c > 0, C„ are suitable constants. 

Note that already the V{tp), cf. (3.3) can be cast in the form (4.5). Furthermore 
the Grassmanian integral P(dV') can be written as the "Lebesgue Grassmanian 
integral" dip'^dil)~ times an exponential 

P(dV) = con5ie--^-/Jo''^-<^-''^?(^«o-|(^»-f^.))V^,- ^^,+^^^- (4 

where dip^dip^ is the Grassmanian integration with "trivial propagators', i.e. with 
the only non vanishing propagator given by {tp^ijj^) = S{(, — 77). Therefore the 
argument of the exponent can also be written in the form (4.5) (involving only 
derivatives of first order at most: as the second derivative can be integrated by 
parts): a property that turns out to be quite important. 

Remark: the difficulty mentioned in the remark following (4.3) shows up very 
clearly here. The integrals in the exponential in (4.5) should be over the rectangle 
[—L/2, L/2] X [0, P] with periodic-antiperiodic boundary conditions. However they 



(4.3) 



(4.5) 



271^3^02000 



are extended to the whole : this means that we have impheitly taken the hmits 
in (3.1). However strictly speaking this does not make sense unless we explain what 
it means to integrate over the whole space a monomial in Grassmanian variables. 
The correct interpretation is the following: 

(a) keeping /3, L finite and using the propagator g{^) in (3.2) one performs the 
integration over the ultraviolet components ij)'^'^^^ and one obtains (4.4) with a 
y(S>) similar to (4.5) but with the coefficient kernels which now depend on L,/?. 

(b) the kernels converge to limits as /3, L ^ oo and the limits Vm°p(^i • . . Cm) 
verify the properties stated after (4.5) (uniformly in (3,L). 

(c) the identity (4.4) means that if one wishes to compute 



/P(dV(>°))^(dV'^-°^)e''('^'V'jf°^">« 
/P((iV(>o))P(#(<o))e^W 

then one can "simply" compute 

/P(dV(^°))e^'<'" 



(4.7) 



(4.8) 



by developing the in powers of the fields and then apply Wick's rule with 

propagator g''-^^ obtaining in this way a combination of integrals of products of 
the kernels in (4.5). And the series converges. 

(d) all divergences due to the infinite extension of the domains of integration 
over the ^ variables disappear when one considers the ratio in (4.3) because of 
the fast decay of the kernels. As one may suspect there will also be exchange of 
limits queries, which are solved again by using the uniformity of the estimates. 
See [BGPS94], for instance, and [BMOOa]. 

Performing the "remaining integral" is therefore similar to the original problem 
except that now the expression V^^^ is more involved but the propagator is simpler 
(being "just" i.e. a propagator with an ultraviolet cut off at scale P(7^). 

In fact the decay of the coefficients Vm^ make the V^'^^ "essentially local" and the 
"remaining integral" is not really harder than the one obtained by replacing y^'^' 

by 



y(o)(^(<o)) = ^„ /%(<o)+^(<o)-^^+ 
J-L/2J0 

/L/2 pl3 



-L/2 . 

where w is a smooth potential with range Pq^ ■ 

In fact in most treatments the above analysis is considered "trivial" and one just 
poses the problem of studying the ratio of integrals 



/P(dV^(^o))e-^""(^'-°')V.g°)"^ 
/P(d?/'(<o))e-^""(V'(s<')) 

with y(°) given by (4.9). 



(4.10) 



It should be noted that the ultraviolet problem does not even arise when the 
fermions are supposed to be located on a lattice because in that case there is a 
minimum length scale, hence a maximum momentum. 

For this reason I concentrate here on the infrared problem, modeled by (4.10). 
However things are not so simple in the case of the Luttinger model: in that model 
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the inverse propagator diverges at large momenta linearly rather than quadrati- 
cally. This generates a non trivial ultraviolet problem: it can be treated in a way 
analogous to the above, but one has to exhibit various cancellations because the 
expressions for the kernels of the effective potential V^^^ are given by apparently 
non convergent integrals. The analysis for the (harder) Luttinger model case is 
carried out in detail in [GS93] . 

§5. The infrared problem and quasi particles. 

The infrared problem (i.e. understanding the large distance properties of the 
Schwinger functions, or equivalently the singularities at finite momentum of their 
Fourier transforms) is more interesting and rich in structure. A naive application 
of perturbation theory leads to facing the fact that the propagators oscillate on 
scale and decay slowly at infinity. 

Quasi particles arise when one attempts to disentangle the oscillations and the 
decay at oo. Technically one remarks that the propagator for the functional inte- 
grals with respect to P{d'ip) can be written, setting ^ = {t, x) and k= {ko, k) in 
general dimension > 1, as 



{27tY+^ -iqo + (£2 - pj,)/2m J 

-J eik,uj) i2nr+^ 

where w is a unit vector and do; is the integration over the miit sphere in R'^ 
normalized to 1 (if the dimension d = 1 integrating over lu means 2^^ '}2iu=± ')• 
An elementary calculation shows that (at least if d is odd) g{x,t, lu), which is 
not unique, can be defined so that e{k, UJ_) — —iko + vpui. ■ k + 0{k^) (see §5 and 
appendix A in [BG90]). 

This means that the free fermion system in a ground state with Fermi momentum 

at pf can be considered as a system of "quasi particles" in the vacmim carrying an 
'Hntrinsid^ linear momentum equal to a Fermi sphere momentum pp w in addition 
to the ^'externaF momentum k. The dispersion relation is almost linear in the 
sense that the system on large scales, i.e. k small, will show a dispersion relation 
essentially identical to e{k, iv) = u ■ k vf- this property seems to remain valid 
even in presence of interaction so that the intuition can be led by the idea that 
the quasi particles must be taken seriously, sec sec. 5 in [BG90]. Formally (5.1) 
can be regarded, c/. also (1.2), as the propagator of a composite field defined as 

i,% = j duie^'^-^^i^%^^^ (5.2) 

with the fields '0x,t,w being Grassmanian fields with propagators g{£,, uj). 

Remark: It is important to stress that the above fields are just Grassmanian fields 
rather than the usual fermionic fields. If we wanted, instead, to think of the quasi 
particles fields as Grassmanian fields which correspond to physical particles it would 
be necessary thai the propagators g(^, uj_) have the "reflection positivity" property, 
[Si74]: and this would also be a criterion to fix the arbitrariness in the choice of 
the representation (5.1) mentioned above. However it is not known whether such 
a choice is possible. The success of the quasi particles-based view of the theory of 
the ground states of fermions leads us to feel that this might be possible. 

It is difficult to give a meaning to the simple and captivating approximation 
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^-i{kot+kx) 



(5.3) 



(27r)'^+i -iko + vfuj_- k 



if > 1 because it corresponds to a system of fermions with a linear dispersion 

relation which therefore will give rise to ultraviolet instabilities: possibly when 
the potential is repulsive and certainly when it is attractive. Note, however, that 
it makes sense if d = 1 because in this case the system is stable even with a 
linear kinetic energy: in fact it becomes almost precisely the Luttingcr model (it 
differs from it because of the presence of "extra" cross terms in the interaction, 
cf. comment following (1.8)). 

However (5.3) gives a good representation of the propagators for sm,a,ll k. There- 
fore the approximation can be expected to be reasonable in the sense that its 
version with an ultraviolet cut-off, i.e. 



can be used to study the infrared problems presented by the (4.9), (4.10). 

The papers [BG90] and [BGPS94] consider the integrals in (4.10) with a propa- 
gator (5.4) and an interaction V^^^ even simpler than the one in (4.9), namely 



where several crossed terms have been eliminated and the interaction has been 
made strictly local (which, in the spinless case, is possible only if we think of the 
system as made with quasi particles because the exclusion principle makes local 
interactions which are quartic in spinless fermion fields vanish identically). 

In the quoted papers it is shown in detail that if d = 1 the understanding of 
the integrals in (4.10) with a propagator (5.4) and an interaction V^^^ given by 
(5.5) suffice to solve in a mathematically complete way the problem of the ground 
state of (1.1) and to obtain Tomonaga's main result that the anomaly at the Fermi 
surface is not and in fact it is an analytic function of the coupling A. Of course 
this is no surprise because of the works [Lu63], [ML65]. 

Why to redo Tom.onaga/s work in m,athem,atically precise way? the point is that 
the notion of "understanding" and of "proof" of a physical result evolve. And many 
of the results that were considered established at the time (~ 1950) have come 
under scrutiny and have been cast into a more rigorous form. This is necessary 
mostly because the attempts to extend them to other cases (namely higher dimen- 
sion or even one dimension with spin, or just with more structure, see question (2) 
in §8) have failed and therefore it becomes necessary to have a clear idea of what 
really fails and what can still be usefully taken over to attack harder problems. 
This necessity is similar to the need that became evident in the 1960's to have 
a more rigorous foundation of the theory of ensembles in equilibrium statistical 
mechanics. It should be clear that revisiting Tomonaga's theory in no way implies 
that there are faults in the original work: it simply did not deal with questions 
that at the time were, rightly, not considered important. 
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(5.4) 




(5.5) 




§6. Renormalization group and the infrared problem. 

We consider the problem of computing the "partition function", i.e. the problem 
of studying the integral in the denominator of (4.9) with a propagator given by 
(5.4) and an interaction given by (5.5). We define 



7^- 1 le^(/cct+u...) ^ 2'^g('\2%u;) 

(ztt)^ —iko + VFtJk 

(6.1) 

h = 0,-1, -2,...,oj = ±,K= (fco, k), = (fcg + vIP)/pI Therefore if Zq 1 



^9^^'\^,u.) = yY.^^ 9^'^ {2% CO) H/ ^5 W (^, c.) + ^5^^-^) ^) (6.2) 

h—O 

and the integration (4.10) can be thought of as an integration over quasi particle 
fields ^l^p which are decomposable as the sum of two fields ip^^^ and 'tp^^^^'' with 
propagators ■^g^°\^,uj) and ^g^^-'^\^,u;): 



I P(#(^o))e^""(v^^) ^ I Pz„(d^W)Pz„(#(^-'')e^'"(^(^'°'+*'-"'» 

(6.3) 

This is convenient because "we know" how to compiite the integral over Vj(") by 
perturbation theory. This means that the techniques to study the integral although 
not trivial are, nevertheless, well established mainly via the results in [Le87] which 
provide us with a technique designed to take advamtagc of the fermionic nature 
of the fields. The result is the expected one: one can basically compute to second 
order of perturbation in A and neglect the rest. More precisely one can prove that 
the result of the integration (in the sense discussed in the remark following (4.5)) 
is 



/ 



Pzo(#(=^-'^)e^""(^*'-"'\ with V^°\y%i^^^-'^) = 

(6.4) 



n>0,a 

where is cither 'ip'^'' , (t^ = ±, or a derivative of this field; and a denotes the labels 
cTj as well as the labels necessary to identify which of the fields $ are differentiated 
and which are not; we also suppose that the tjj'^ fields are to the left of the V~ 
(recall, however, that the fields are Grassmanian so that they anticommute) . The 
number n must be, obviously, even. 

The kernels Va{^i, . ■ . ,^n) are analytic functions of Ao, I'd convergent with a con- 
vergence radius which is independent of the sizes L,I3 of the system and decay 
exponentially fast on scale p^^, in the sense following (4.5): the technique for this 
proof is in [Le87] and the analysis can be found in [BGPS94]. 

It would not be wise to simply iterate the procedure calling V^~^'' the "effective 

interaction'' y*""' and (with the procedure already used in the previous splitting, 
of. (6.2)) splitting the field ip^-~^^ into 'tp^~^^ +%p^-~^^: because the result, under 
further iteration, would be a progressive worsening of the bounds on the constants 
and a decrease in the convergence radius of the expansions in powers oi]lQ,UQ. This 
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is basically due to the fact that i/q cannot be arbitrary under the only c;ondition 
that it is small: it has to be tuned so that the Schwinger functions have singularities 
at momentum ±pf- 

The correct procedure will be to adjust rcnnu'sively the value of as a function of 
Ao keeping the singularity at fixed momentum. For this purpose one distinguishes 

the relevant part of the interaction V^^^ from the irrelevant part. The latter is not 
small (in a sense it is the most important part!): the relevant part is just a part 

of V^^"' whose absence would allow us to perform a naive recursion without loss in 
the size of the constants or of the convergence radii. 

The identification of the relevant terms is a matter of dimensional analysis, at 

least in the simplest cases, sec [Ga85], [BG95]. What we call here "relevant" is a set 
of terms (see below) which in the usual nomenclature of the renormalization group 
approaches is further divided into relevant and marginal terms. In the present case 
the relevant part is the "/ocaZ partr of the terms which are quadratic and quartic 
in the fields. If da,b denotes a Kronecker delta such local part is defined by 



j=i iD-o; 

where = [di, 9,;). 

Applying the operator L to the expression in (6.4) one realizes that the result 
can be espressed as a linear combination of the following Grassmanian monomials 

P2=- [ E<«C,«'^^. (6-6) 

d^ 



Note that there is only one possible non zero local term of fourth order in the field 
because of the Fermi statistics and the facts that w = ±1 and that our fermions 
are spinless. Note also that the above localization operation would be completely 
diff'erent {and useless) if we had not introduced the quasi particles: for instance Fx 
would simply vanish because of the Fermi statistics as it would involve the square 
of Fermi fields. 

Applying the operator £ to v'^^\\/Zoil)'^-~^^) one obtains an expression 

A^°^ Zo^Fi(^(^-i)) +fWZoF2(V(^-^))+ 

+ a(o)ZoF3(V^--^^) + C^°^^oi^4(V(--^^) + (1 - >C)F^°^ 

where A , C, are simple combinations of integrals of the kernels, see 

(6.4), for the monomials quadratic or quartic in the fields and are therefore analytic 
in Xo,vo. 

The simplest way to procecid is to consider a representation of Pz„ (dV') analo- 
gous to (4.6) for the integrations over the quasi particle fields, with the quadratic 
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form i>f-'^+{dt + i-^dl - ^p%))4^-^^- replaced by Ea,=± V-g'^^^l^t - 
ivpLodx)lL^-i{d)iplj^^ where Th{d) is the operator which multiphes the Fourier 
transforms by Tii{)i) = e+(^ (k„+vpk ). Jj^jJqqjJ choice attributes to the 
field the correct propagator ^'-"^'(O) c/. (6.2). Then the integral to be 

performed looks like 

const J e'^" ^-=± / ^-9)r_i(a 



e 



■ e 



dtp~^dtp 



and it is natural to collect the quadratic parts defining 

Z_i Zq + aoZo, C-iZ-i (C^"^ - ao)-^o 



(6.9) 



Here the factor 2 in front of the "relevant coupling" Vqis natural: if the theory is 
developed only to first order (a deed called "power counting" in the jargon of the 
renormalization group, see eq. (6.18) in [BG95]) the result is that Z_i = Zq and 
= 21/0. This transforms the integral into 



/ 



(6.10) 



The procedure can now be iterated and the integrals can be performed in the 
same way defining recursively Xh,i/h,Ch called "running couplings", Zh called 
"wave function renormalization" and Wh called "irrelevant operator^ for h = 
0,-1,-2 The ratio Zh/Zh+i will be called the "wave function renormaliza- 
tion rate". 

Remarks: (i) Note that 

Whi^/zTrii^^^-^))) = (1 - £)F^°^(yZoV(^-^^)+ 



(6.11) 



(ii) In [BG90], [BGPS94], [BG95] the definitions of Wh and/or of the localization 
operator C are slightly different: more involved but, possibly, more convenient for 
performing estimates. The choices differ, however, by "irrelevant terms" and are 
(therefore) equivalent. 

The basic bound is that the running couplings on scale h and the kernels that are 

coefficients in the definition of Wh are analytic in V/^, = (Xh' , i^h' , Ch') for h' > h 
under the condition that all the latter quantities and \Zh'/Zh'i — 1| for h' > h are 
small enough independently on the size L of the system and also on the scale h. 

The relation between the running couplings can be written, therefore, as 
Vh-1 = Mh{vh + Bhivh^Vh+i,...,Vo;—^,...) , 

Zh - Zh ^'-'^^ 

1= 1^ + Ahivh, Vh+i^ ■ ■ ■ ^ Ho'^ ■ ■ ■)) ' 

^h-l 
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where Mh is a diagonal matrix with matrix elements {Zh/ Zfi^i)^ , 2 Zh/ Zh-i, 
Zh/Zh-x and the functions Bh,Ah are analytic under the conditions that the 
running couplings and the renormalization rate are small enough independently of 
h < 0. 

Furthermore one can prove, [BG90], [BGPS94], that after eliminating the wave 
function renormalization rates in the first of (6.12) by using recursively the second 
equation the relation (6.12) can be written 

^^=l + B',iv,^,,v„...,v,) ^^^^^ 

llh = ^ Vh+i + ^ft(2^ft+i, Vh,---, Ho) 

with i?^, B_i^ analytic in their arguments u^/, under the conditions that the run- 
ning couplings are small enough independently of /i < 0, and A is the diagonal 
matrix with diagonal (1,1,2). 

At this point the strategy is clear: the (6.13), called the beta functional, defines 

a "flow" , called the renormalization flow, in the space of the running couplings 
V = (A, v, Q which is well defined as long as the running couplings themselves 
stay small enough in the above sense so that perturbation theory can be applied 
to pass from one scale to the next {i.e. as long as 1 1; /J is so small to be inside the 
convergence domain of the series expressing the functions 5^, B_f^). The initial 
values are Aq, i^O) Co: Aq will have to be taken small (so that the perturbation theory 
can be applied at least to perform the first integration, i.e. the integration over 
tp^^^), also will have to be taken small and Co must be taken Co = (because we 
start the analysis from (5.5) which contains no term like F^, see (6.6.). 
Note that if we studied the problem without the approximation introduced after 
(5.5) (starting from the Hamiltonian (3.3)) we would end up after the ultraviolet 
integration with a V^^^ which is much more involved than the expression (5.5) with 
which we began the infrared analysis. Therefore the localization operator acting 
on y*^") may (as indeed it does) produce also a term of the type i^4: which implies 
that Co will have a non zero value. This however does not help to provide us with 
more freedom because this Co is not arbitrary being a function of the couplings A, v 
in the original Hamiltonian (3.3): it is intrinsic to the model that there are only 
two adjustable constants in the Hamiltonian. 

§7. The vanishing of the beta function: the role of the Luttinger model. 

The only freedom that we have to "make things work" is the choice Aq (or of A 
in the theory without approximations) small "enough" and the selection of uq (or 

of V in the theory without approximations). 

What has been sketched until now is technically involved, but it is not really diffi- 
cult because the techniques for a rigorous multiscale analysis has been established 

since a long time. In many problems, among which the rf-dimensional Fermi sys- 
tems, see [BG90], or the scalar field theories, see [Ga85], [BG95], it is, in a sense, 
a pure matter of technical routine work to prove that the successive integration 
of the components of the fields on various scales can be reduced to the study of a 
renormalization flow like (6.12), under the assumption that the running constants 
stay small enough. 

In itself such a work does not really mean much unless one is able to exhibit a 
trajectory of the flow that keeps the running constants small enough to allow us 
to apply perturbation theory to compute the effective potentials from one scale to 
the next via convergent series for the kernels coefficients. 

Remarks: (1) It is iii(lec;(l very easy to work out alternative, substantially different, 
integration algorithms which work as beautifully as the one described above but 
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which are completely useless: a nicx' and relevant example is obtained by keeping 
Zh = 1 and, correspondingly, not eliminating the quadratic terms with coefficients 
ah that are generated at each integration of a field component i/;^''^ Everything 
works fine: however one can prove that no initial data Ao,z/o exist which keep the 
constants Xh, fhi ochi Ch small for all h if ao = 0,(^o = (or, in the theory without 
approximations, if A and v are the only free parameters )! 

(2) The only interest of the latter remark (1) lies in the fact that if it had been 
possible to find vq as a function of Aq so that Xh,i^h,cth, Zh stayed small for all < 
then it would follow (with some extra work, see [BGPS94]) that the singularity of 
the one particle Schwinger function at the Fermi surface would be a disc;ontinuity, 
as in the non interacting case. This failure shows that the latter property cannot 
be derived as a consequence of a perturbation analysis, if true (it is not! as we 
know since Tomonaga's work). 

(3) Of course the impossibility of a normal Fermi surface in generic 1-dimensional 
many fermions systems has been known since the early days of many body theory 
as a consequence of the divergence of the second order eorrcictions to the Schwinger 
functions (divergence of the "self-energy" which occurs in dimension 1). 

Therefore after the above (long and tiring) set up of the renormalization flow the 
real work starts. 

It was well known that there were interesting conjectures about a number of 
cancellations that occur in the theory of one dimensional spinless fermions. In our 
language they can be summarized by saying that "the beta function effectively 
vanishes", see [So79]. 

It is very difficult however to find a proof by computing the beta functional to all 
orders. Introducing the function, called the "beta function" , 

Hiv)"^ \vca B^{v,v,...,v), v={X,u,OgR' (7.1) 

— n— ^— oo 

a major simplification occurs by realizing, see [BG90], [BGPS94], that there is a 

solution to the flow generated by the recursion (6.13) with | so small that all 
running couplings stay within the convergence domain of the functions pro- 
vided the components (relative to the contants X and /3i,/32 of the beta function 
^ = (/?i,/32,/?3) vanish. 

A proof of the latter property looks like an easier problem. However the first 
proof came from an indirect argument that we describe below. 

(1) one first repeats the above analysis for the Luttinger model: this was assumed 
possible in [BG90], [BGM92], and later proved in [GS93]: the difficulty being 
essentially in the treatment of the ultraviolet problem. 

(2) once the beta functional has been defined for the Luttinger model it has been 
remarked in [BG90] that although the beta functionals of the Luttinger model 
and the corresponding one of the "realistic" model in (1.1) are different the beta 
functions of the two models coincide, [BG90]. 

(3) one checks that if 3i or /32 did not vanish this would contradict the exact 
solution of the Luttinger model, [BG90], [BGM92], [BGPF94]. 

This allowed to prove the existence of a trajectory of the renormalization flow 
which stays close; to the origin and within the radius of convergence of all the 
expansions in the running couplings considered (including that of the beta func- 
tional) thus completing the theory of the 1-dimensional spinless Fermni gas at 
small coupling yielding also the analytieity of the anomaly exponent in A. 

The third component of the beta function does not vanish: this is however not 
necessary. In fact the constant Uh tends to diverge and very fast (0(2~'')): but 
this can be easily counterbalanced because the initial value of i/q is free and we can 
tune it so that Vh does not diverge. In fact there can be only one value of which 
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has this property (essentially beeaiise, near a hyperbolic equilibrium attractive or 
marginal along the A-axis and repulsive along the i^-axis, given A small enough 
one can find only one value of v small such that the evolution of the datum A, v 
stays close to the origin forever). 

§8. Outlook. 

One can wonder whether a simpler symmetry argument can be given to prove that 
the beta function /3i,/32 components vanish: after [BG90] this has been attempted 
in [MD93] who however neglect the (necessary) presence of the decreasing cut-off 
I^Vq, for h — > — oo, arising in the successive integrations of the infrared scales. 

Nevertheless such a proof must be possible and recently new attempts at finding it 
are being considered, [BMOOb] : this would be a major achievement which, however, 
would not diminish the importance of the idea of Luttinger to enucleate out of 
Tomonaga's model a much simpler model that rc;ally catches all its main features. 
The existence of the indirect proof of the vanishing of the beta function based in 
an essential way on the [ML65] solution of the Luttinger model, is an important 
factor in the search of the proof based on symmetry considerations. 

Mainly through the work of Mastropietro and collaborators the Luttinger model 

(together with its extension by Mattis, [Ma64]) has received several other appli- 
cations in the frame of the renormalization group approach to the ground state 
of Fermi systems; for a complete review see [GMOO] where, among other results, 
the interesting phenomenon of the A dependence of the anomalies in the higher 
Schwinger functions is summarized as well as the results of the theory of one dimen- 
sional fermions in a periodic or quasi "periodic potential (with Fermi momentum 
PiT in the bands or at the top of a band). 

Another somewhat unexpected development has been the determination of the 
asymptotic behavior of the correlation functions in eight vertex models or XYZ 
models in a magnetic field for values of the parameters that do not correspond to 
exactly soluble points. [Ma99b], [BMOOa]. A corollary of the latter papers should 
be the determination of the critical exponents of Ising type models even with 
non nearest neighbour interactions, see [Ma99b]. Applying the remark in [Sp99] 
concerning the translation of the next neighbor 2 dimensional Ising model into 
a fermionic spinless 1-dimensional problem it could lead to a strong extension 
of the theory of the critical point in such model covering a wide range of other 
models which (unlike the case considered in [Sp99]) show non universal critical 
exponents, see also remarks to sec. 1.6 of [BMOOa], [MaOO] and the review [GMOO]. 
Furthermore the technique used in [Ma99b], [BMOOa] is quite different from that 
in [Sp99]. 

I conclude by mentioning a few open problems for small coupling: 

(1) spinning fermions on a line with attractive interaction. See [BM95] for the 

repulsive case 

(2) spinless fermions on two parallel lines with attractive interaction, see [Ma99a] 
for the repulsive case. 

(3) two and three dimensional Fermi surface properties(!). 

(4) a mathematical proof of the identity p = irpp between the density p and the 
Fermi momentum p/^- , see [Lu60] and §3 above, at small coupling A. 

and it is, perhaps, surprising to find the first two in a list of open problems, 
while (3) is well known to require new ideas to be understoon and (4) seems 
understandable with the techniques that we know. 

We can say that the work of Luttinger [Lu63] proposing the exactly soluble model 
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bearing his name and providing some of the ideas that contributed to its exact so- 
lution in [ML65] has been, and remains, a landmark in Condensed Matter Physics 
(with its clarifying function of Tomonaga's work) and in Mathematical Physics 
(as an example of how to build a model that catches all the relevant features of 
a realistic model and, yet, it is more tractable) and it has shed new light on the 
difficult subject of the theory of low temperature quantum systems. 
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